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ABSTRACT 



We discuss the forces on grains exposed to anisotropic 
radiation fields, including the usual "radiation pres- 
sure" force and also recoil forces due to photoemit- 
ted electrons, photodesorbed hydrogen atoms, and 
hydrogen molecules which form on the grain surface. 
We show that these forces can lead to grain dynamics 
in photodissociation regions which result in enhanced 
dust-to-gas ratios. Since the gas heating is probably 
dominated by photoelectric emission from dust, this 
might explain the unusually high gas temperatures 
inferred from ISO observations of molecular hydro- 
gen in photodissociation regions. 
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1. INTRODUCTION 



be a major mechanism, but calculated heating rates 
are inadequate. Here we consider the possibility of 
enhanced dust-to-gas ratios in PDRs, which would 
imply increased photoelectric heating rates. 



2. GRAIN DYNAMICS 



Enhanced dust-to-gas ratios can result when grains 
drift with respect to the gas. We will begin by de- 
scribing the motion of the gas (see Figure 1). Ion- 
izing radiation from the hot star located near the 
molecular cloud (MC) photoevaporates material at 
the cloud surface, which then expands into the in- 
tercloud medium in a photoevaporative flow (PF). 
Thus, the ionization front propagates into the cloud. 
In our analysis, we adopt the frame of reference in 
which the ionization front is stationary; in this frame 
material flows through the PDR. 



Far ultraviolet radiation incident on a cloud of molec- 
ular material dissociates molecules, giving rise to a 
surface layer of largely atomic material, called a pho- 
todissociation region (PDR). When a hot star is lo- 
cated near a molecular cloud, much of the star's ra- 
diant output is absorbed and reprocessed in the re- 
sulting PDR. Dust in the PDR efficiently absorbs UV 
photons, which results in thermal IR emission. The 
warm gas cools via forbidden transitions in metal 
atoms and ions and via rovibrational transitions in 
H2. Thus, the conditions in PDRs are probed ob- 
servationally by IR studies. Many theoretical PDR 
models have been developed, and include quantita- 
tive determinations of IR diagnostics (e.g. Tielens 
& Hollenbach 1985; Draine & Bertoldi 1996; see also 
the review by Hollenbach & Tielens 1997 and refer- 
ences therein). The detailed understanding of PDRs 
will yield information on the structure of molecular 
clouds and their interaction with young stars, and 
thus ultimately on the process of star formation. 

As discussed by Draine & Bertoldi in this volume, ob- 
servations with the ISO short wavelength spectrome- 
ter have yielded the H2 rotational distribution func- 
tions in PDRs, from which surprisingly high values of 
gas temperature (500 - 1000 K) have been inferred. 
It is not yet clear what heating processes are able to 
maintain the gas at such high temperatures. Pho- 
toelectric emission from dust grains is expected to 




Figure 1. The physical picture. 

If no other forces acted on the grains, they would sim- 
ply be dragged along by the gas. However, there are 
forces associated with the anisotropic stellar UV ra- 
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diation field. In addition to the usual "radiation pres- 
sure" force due to absorption and scattering (-F ra d), 
there are also recoil forces, due to photoemitted elec- 
trons (F pe ) and due to photodesorbed adatoms and 
H2 molecules which form on the grain surface (-Fpd)- 
As a consequence of these "radiative" forces, the 
grains drift with respect to the gas, at speed thrift • If 
the gas moves at speed v with respect to the ioniza- 
tion front then the grains move at speed vq — Wdrift- 
Deep within the PDR, where the stellar UV has been 
attenuated, thrift = 0. The grain overdensity near 
the ionization front is obtained using the continuity 
equation: 

= 11- ! , (1) 



n 



vo 



where no (rii) is the grain number density deep within 
the PDR (near the ionization front) and thrift, i 1S the 
drift speed near the ionization front; we assume one- 
dimensional flow and ignore magnetic fields. Since 
the drag force is roughly proportional to the drift 
speed, the overdensity is also given by 



n 



stop 



(2) 



where F tot is the sum of the radiative forces near 
the ionization front and F stop is the magnitude of 
the radiative force which would result in thrift = v o 
(in which case the grains would be stopped near the 
ionization front). 



3. EVALUATION OF RADIATIVE FORCES 



ionization potential depends on the grain charge; the 
charge distribution is determined by the rates of elec- 
tron attachment by accretion from the gas and re- 
moval by photoemission. Thus, the average photo- 
electric force depends on the ambient conditions; in- 
cluding the spectrum and energy density of the radia- 
tion and the gas density, ionization, and temperature; 
through their effects on the charge distribution. For 
a fixed charge state, the force is given by 



F pc = 7ra 2 / dv —^AYQ ahs S\/2m e E eh 



(4) 



where v^ — IP/h, Y is the photoelectric yield, the 
emission asymmetry factor A(hv, a) measures the 
asymmetry in the emission of photoclcctrons over the 
grain surface, the recoil suppression factor S accounts 
for electron emission in directions other than the sur- 
face normal, and m e is the electron mass. 

We follow the simple prescription of Kerker & Wang 
(1982) for determining the asymmetry factor A. The 
probability of photoemission from any site on the sur- 
face is taken to be proportional to the electric inten- 

sity |E| just below the surface at that point. Thus, 



A(hv, a) = 



X sin 9 cos 9\E{9)\ 2 d9 
f;sinfl|E(0)| 2 ^ ' 



(5) 



where 9 is the polar angle with respect to the direc- 
tion of the incident radiation. In Figure 2 we display 
A(hv, a) as a function of incident photon energy for 
graphite grains of various sizes. 



Throughout, we assume spherical grains, so that Mie 
theory can be applied for the grain optics. 



3.1. Radiation Pressure Force 

The radiation pressure force is given by 

F rad = no 2 / dv u v [Q a bs + Q SC a(l - (cos 9))\ , 
Jo 

(3) 

where a is the grain radius, vn = 13.6 eV/h, u v is 
the radiation energy density per frequency interval, 
and (cos 9) is the usual scattering asymmetry factor. 
We compute the absorption and scattering efficiency 
factors, Qabs and <5 sca respectively, using a Mie the- 
ory code derived from BHMIE (Bohrcn & Huffman 
1983) with dielectric functions as described by Draine 
& Lee (1984) and Laor & Draine (1993). 




hv (eV) 



3.2. Photoelectric Force 



We adopt a physical model for the photoelectric emis- 
sion process (see Bakes & Tielens 1994, Weingartner 
& Draine 1999a). Of course, the photoemission rate 
and the total energy of the emitted electrons depend 
on the grain's ionization potential, IP; for simplic- 
ity we take the electron energy E c \ = (hv — IP)/2, 
where hv is the energy of the absorbed photon. The 



Figure 2. The emission asymmetry factor A as 
a function of the incident photon energy hv, for 
graphite grains. The grain radius is indicated for each 
curve. 

In determining the recoil suppression factor S we as- 
sume that the electrons emerge symmetrically with 
respect to the local surface normal, with a "cosine- 
law" angular distribution (i.e. the emission rate 
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at angle ip with respect to the surface normal oc 
sin xjj cos ip). For an uncharged grain, this would im- 
ply S = 2/3. When the grain is charged, the electron 
escapes on a hyperbolic trajectory, so that when it is 
at infinity its velocity vector makes an angle with 
respect to the surface normal which differs from the 
corresponding angle at the surface, ip. See Weingart- 
ner & Draine (1999b) for the detailed calculation. 



3.3. Photodesorption Force 

In our simple model, we consider the exchange of 
hydrogen between the gas phase and the grain sur- 
face. We assume that the surface is entirely covered 
with chcmisorption sites, as would be the case for 
bare graphite or silicate grains without ice mantles. 
When a gas phase H atom strikes the grain surface, 
it arrives either at an empty site or at a site which 
is already occupied by an adsorbed H atom. In the 
former case, it occupies the site with some sticking 
probability; otherwise it reflects and remains in the 
gas phase. In the latter case, it recombines with the 
resident H with some recombination probability; oth- 
erwise it reflects and remains in the gas phase. We 
assume that H2 formed on the grain surface is im- 
mediately ejected from the grain. Adatoms are re- 
moved from the grain surface via recombination and 
photodesorption. 

Suppose that the grain's spin axis is aligned with the 
direction of the radiation anisotropy. Since the pho- 
todesorption rate is greater in the illuminated hemi- 
sphere than in the non-illuminated hemisphere, rela- 
tively more H leaves the surface as atoms in the for- 
mer and as molecules in the latter. A net recoil force 
is expected, for two reasons. First, the total number 
of particles leaving the illuminated hemisphere ex- 
ceeds that for the non-illuminated hemisphere. Sec- 
ond, the outgoing atoms likely carry a different mo- 
mentum than the outgoing molecules. A net recoil 
force of zero is possible, but only if the momenta and 
sticking and recombination probabilities conspire to 
have the right values. The photodesorption force is 
greater when the spin axis makes a non-zero angle 
with respect to the radiation direction, since H which 
arrives at non-illuminated sites can be transported to 
illumination and photodesorbed. For a fuller discus- 
sion of the adopted values for relevant parameters 
and the details of the calculation, sec Weingartner & 
Draine (1999b). 
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Figure 3. Ratios of the radiative forces to F st0 p, the 
force needed to stop graphite grains relative to the 
ionization front, near the ionization front. 
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4. THE S140 PDR 



Timmermann et al. (1996) studied the PDR asso- 
ciated with the S140 HII region using ISO and con- 
cluded that the gas temperature is £ 500 K. We have 
evaluated the radiative forces for graphite and silicate 
grains and the conditions in this PDR. To determine 
-Fstop , we employed the expression for -Fdrag in Draine 
& Salpeter (1979) and evaluated vo using a highly 
simplified model for the flow in the PDR; see Wein- 
gartner & Draine (1999b) for details. We display the 
ratios of the radiative forces to F stop , as a function 
of grain size, in Figures 3 and 4. 



Figure 4. Same as Figure 3, but for silicate grains. 



5. CONCLUSIONS AND FUTURE WORK 



For the S140 PDR, the forces associated with the 
anisotropic radiation are apparently large enough to 
result in significant grain overdensities. It is impor- 
tant to note, however, that the total radiative force 
is highly uncertain, due to uncertainties in the mi- 
crophysics of the photoemission and photodesorption 
processes. 

Our program of study will include an examination 
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of more PDRs and also simulations of the dynam- 
ics which result when F tot > F stop ] i-e. when the 
grains near the ionization front are actually pushed 
deeper into the PDR. In this case, the further up- 
stream a grain is pushed, the less effective the push, 
because the radiation is attenuated by other grains 
downstream. Also, we will add magnetic fields to 
the analysis. Finally, of course, we will estimate the 
implications for gas heating. 
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